First Results from the XENON10 Dark Matter Experiment 
at the Gran Sasso National Laboratory 
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The XENON10 experiment at the Gran Sasso National Laboratory uses a 15 kg xenon dual phase 
time projection chamber (XeTPC) to search for dark matter weakly interacting massive particles 
(WIMPs). The detector measures simultaneously the scintillation and the ionization produced by 
radiation in pure liquid xenon, to discriminate signal from background down to 4.5 keV nuclear 
recoil energy. A blind analysis of 58.6 live days of data, acquired between October 6, 2006 and 
February 14, 2007, and using a fiducial mass of 5.4 kg, excludes previously unexplored parameter 
space, setting a new 90% C.L. upper limit for the WIMP-nucleon spin-independent cross-section 
of 8.8 x 10 _44 cm 2 for a WIMP mass of 100 GeV/c 2 , and 4.5 x 10" 
30 GeV/c 2 . This result further constrains predictions of supersymmetric models. 

PACS numbers: 95.35.+d, 29.40.Mc, 95.55.Vj 



The well-establisbed evidence for non-baryonic dark 
matter [l|, 0, 9 is a striking motivation for physics be- 
yond the Standard Model of particle physics. Weakly 
interacting massive particles (WIMPs) [J] as dark mat- 
ter candidates arise naturally in various theories, such 
as Supersymmetry, Extra Dimensions, and Little Higgs 
models @ H, 0, S- Since by hypothesis the WIMPs in- 
teract through the weak interaction and can efficiently 
transfer kinetic energy by elastically scattering from 
atomic nuclei, the WIMP model can be tested by search- 
ing for nuclear recoils in a sensitive, low-radioactivity 
detector [H, [nj. Predicted event rates are less than 
0.1 events/kg/day, with energy depositions of the order 
of 10 keV. 

XENON10 is the first 3-D position sensitive TPC de- 
veloped within the XENON program to search for dark 
matter WIMPs in liquid xenon (LXe) [ll|. Dual phase 
operation enables the simultaneous measurement of di- 
rect scintillation in the liquid and of ionization, via pro- 
portional scintillation in the gas [T3 |. The ratio of the 
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two signals is different for nuclear (from WIMPs and neu- 
trons) and electron (from gamma and beta background) 
recoil events [13j . providing event-by-event discrimina- 
tion down to a few keV nuclear recoil energy [3, EH • In 
March 2006 the detector was deployed underground at 
the Gran Sasso National Laboratory (LNGS) [16( , where 
it has been in continuous operation for a period of about 
10 months, with excellent stability and performance [171 ]. 
The TPC active volume is defined by a Teflon cylinder 
of 20 cm inner diameter and 15 cm height. Teflon is used 
as an effective UV light reflector [18j and electrical in- 
sulator. Four stainless steel (SS) mesh electrodes, two 
in the liquid and two in the gas, with appropriate bias 
voltages, define the electric fields to drift ionization elec- 
trons in the liquid, extract them from the liquid surface 
and accelerate them in the gas gap. For the dark mat- 
ter search reported here, the drift field in the liquid was 
0.73 kV/cm. 

Two arrays of 2.5 cm square, compact metal-channel 
photomultiplier tubes (PMTs) (Hamamatsu R8520-06- 
Al) detect both the direct (51) and proportional (52) 
scintillation light. The bottom array of 41 PMTs is in the 
liquid, 1.5 cm below the cathode mesh, to efficiently col- 
lect the majority of the direct light which is preferentially 
reflected downwards at the liquid-gas interface. The top 



2 



array of 48 PMTs, in the gas, detects the majority of the 
proportional scintillation light. From the distribution of 
the PMT hits on the top array, the event location in XY 
can be reconstructed with a position resolution of a few 
millimeters. The third coordinate is inferred from the 
electron drift time measured across 15 cm of LXe, with 
better than 1 millimeter resolution. The PMTs are dig- 
itized at 105 MHz with the trigger provided by the 52 
sum of 34 center PMTs of the top array. 

The TPC is enclosed in a SS vessel, insulated by a vac- 
uum cryostat, also made of SS. Reliable and stable cryo- 
genics is provided by a pulse tube refrigerator (PTR) [19| 
with sufficient cooling power to liquefy the Xe gas and 
maintain the liquid temperature at -93 °C. The Xe gas 
used for the XENON10 experiment was commercially 
procured with a guaranteed Kr level below 10 part per 
billion (ppb). The XENON10 detector is surrounded by a 
shield made of 20 cm-thick polyethylene and 20 cm-thick 
lead, to reduce background from external neutrons and 
gamma-rays. At the Gran Sasso depth of 3100 meters 
water equivalent, the surface muon flux is reduced by a 
factor of 10 6 , such that a muon veto was not necessary 
for the sensitivity reach of XENON10. 

We report here the analysis of 58.6 live-days of WIMP- 
search data taken with the XENON10 detector at Gran 
Sasso during the period between October 6, 2006 and 
February 14, 2007. The analysis was performed "blind'' , 
i.e. the events in and near the signal region were not an- 
alyzed until the final signal acceptance window and event 
cuts were tested and defined, using low energy electron 
and nuclear recoils from calibration data, as well as 40 
live-days of "unmasked" WIMP-search data. In this let- 
ter we interpret the data in terms of spin-independent 
WIMP-nucleon scattering cross-section. Another letter 
will focus on the spin-dependent interpretation of the 
same data [20|. More details on the analysis and the es- 
timation of backgrounds will be reported elsewhere [l7T | . 

A total of ~ 10 electron recoil events, in the a priori 
set energy range of interest (4.5 to 26.9 keV nuclear recoil 
equivalent energy) for the WIMP search, were collected 
with 137 Cs source. The number of events in the final 
fiducial volume is about 2400, 1.3 times the statistics 
of the WIMP-search data. The detector's response to 
nuclear recoils was obtained from 12 hours irradiation 
in-situ, using a 200 n/s AmBe source. 

The 51 signal associated with each triggered event is 
searched for in the off-line analysis. By requiring a co- 
incident signal in at least two PMTs, the efficiency of 
the 51 signal search algorithm is larger than 99%, above 
a threshold of 4.4 photoelectrons (pe), or 4.5 keV nu- 
clear recoil equivalent energy. The 52 hardware trigger 
threshold is 100 pe, corresponding to about 4 electrons 
extracted from the liquid, which is the expected charge 
from an event with less than 1 keV nuclear recoil equiva- 
lent energy [l4[ • The 52 trigger efficiency, with a software 
threshold of 300 pe, is more than 99% for 4.5 keV nuclear 
recoils. Basic-quality cuts, tuned on calibration data, are 
used to remove uninteresting events (e.g. multiple scat- 



ter and missing 52 events), with a cut acceptance for 
single-scatter events close to 99%. 

Energy calibration was obtained with an external 57 Co 
gamma ray source and with gamma rays from metastable 
Xe isotopes produced by neutron activation of a 450 g 
Xe sample, introduced into the detector after the WIMP 
search data taking. The 51 and 52 response from the 
131m Xe 164 keV gamma rays, which interact uniformly 
within the detector, were used to correct ±20% variations 
of the signals due to the position dependence of the light 
collection efficiency. 51 for 122 keV gamma rays, after 
position-dependent corrections, gives a volume-averaged 
light yield, L y , of 3.0 ± 0.1(syst)±0.1(stat) pe/keVee 
(keVee is the unit for electron-recoil equivalent energy) 
at the drift field of 0.73 kV/cm. The nuclear recoil 
equivalent energy (in unit of keV) can be calculated as 
E nr = Sl/Ly/Ceff • S e /S n . Here C e ff is the nuclear- 
recoil scintillation efficiency relative to that of 122 keV 
gamma rays in LXe at zero drift field. We used a con- 
stant C e ff value of 0.19 which is a simple ass ump tion 
consistent with the most recent measurements [2l|, l22j . 
5 e and 5„ are the scintillation quenching factors due to 
the electric field, for electron and nuclear recoils, respec- 
tively. S e and S n were measured to be 0.54 and 0.93, 
respectively, at a drift field of 0.73 kV/cm [Tl ]. 

Background rejection is based on the ioniza- 
tion/scintillation (52/51) ratio, which is different for nu- 
clear and electron recoils in LXe. Figure [T] shows the 
energy dependence of the logarithm of this ratio for elec- 
tron recoils from 137 Cs gamma ray and for nuclear recoils 
from AmBe fast neutron calibrations. The separation 
of the mean Logio (52/51) values between electron and 
nuclear recoils increases at lower energy. In addition, 
the width of the electron recoil band is also smaller at 
lower energy. The combination of these two effects gives 
a better electron recoil rejection efficiency at the lower 
energy, reaching 99.9%. The different ionization density 
and track structure of low energy electrons and Xe ions in 
LXe result in different recombination rate and associated 
fluctuations, which might explain the observed behavior. 

The Login (52/51) values increase with decreasing en- 
ergy for both electron and nuclear recoils, within the en- 
ergy window of interest (4.5 - 26.9 keV nuclear recoil 
equivalent energy). The same behavior has been ob- 
served previously in small prototype detectors [3, EH- 
In our analysis, we subtract the energy-dependent mean 
Logio (52/51) from the electron recoil band to obtain 
ALogio(52/51) for all events. After this band flatten- 
ing, the energy window of interest for the WIMP search 
is divided into seven individual energy bins (see Ta- 
ble |T|. For each energy bin, the nuclear recoil acceptance 
window is defined to be between ALogio (52/51) = fi 
and ALogio (52/51) = A* — 3cr. Here fj, and a are the 
mean and sigma from a Gaussian fit of the nuclear recoil 
ALogio(52/51) distribution. The nuclear recoil accep- 
tance efficiency is the fraction of nuclear recoil events 
within the acceptance window. The ALogio(52/51) dis- 
tribution for electron recoils from the 137 Cs data is found 
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Nuclear Recoil Equivalent Energy (keV) 

FIG. 1: Logio(S2/Sl) as a function of energy for electron re- 
coils (top) and nuclear recoils (bottom) from calibration data. 
The colored lines are the mean Logio(S2/51) values of the 
electron recoil (upper, red) and nuclear recoil (lower, blue) 
bands. The region between the two vertical dashed lines is 
the energy window (4.5 - 26.9 keV nuclear recoil equivalent 
energy) chosen for the WIMP search. An S2 software thresh- 
old of 300 pe is also imposed (black lines). 



13 Cs calibration data and unmasked WIMP search data. 
These events were identified to be multiple-scatter events 
with one scatter in the non-active LXe mostly below the 
cathode and a second scatter in the active LXe volume. 
The 52 signal from this type of event is from the interac- 
tion in the active volume only, while the SI signal is the 
sum of the two ST's in both the active and non- active 
volume. The result is a smaller S2/S1 value compared 
to that for a single-scatter event, making some of these 
events appear in the WIMP-search window. Two types 
of cuts, one using the SI signal asymmetry between the 
top and bottom PMT arrays and the other using the ST 

hit pattern, defined as SIrms = \J \ 5Z(Sli ~ S^) 2 (* = 
1, n), on either the bottom or the top PMT array, are de- 
fined to remove these anomalous events. The SI signal 
from the scatter outside the active volume tends to be 
clustered on a few of the bottom PMTs (larger SIrms), 
while the SI signal from a normal event in the active vol- 
ume is distributed more evenly over the PMTs (smaller 
S1_rms)- A large fraction of events that leaked into 
the WIMP-signal window are of this type of background 
and could be removed by the cuts discussed above. The 
cut acceptance e c for single-scatter nuclear recoil events, 
based on AmBe neutron calibration data, is listed in Ta- 
ble m 



TABLE I: The software cut acceptance of nuclear recoils e c , 
the nuclear recoil acceptance A nr , and the electron recoil re- 
jection efficiency R er for each of the seven energy bins (E nr 
in nuclear recoil equivalent energy). The expected number 
of leakage events, Ni ea k, is based on R er and the number of 
detected events, N ev t, in each energy bin, for the 58.6 live- 
days WIMP-search data, with 5.4 kg fiducial. Errors are the 
statistical uncertainty from the Gaussian fits on the electron 
recoil ALogio(S f 2/5'l) distribution. 
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empirically to be statistically consistent with Gaussian 
fits, except for a small number of "anomalous leakage 
events" . From these fits, we estimate the electron recoil 
rejection efficiency and predict the number of statistical 
leakage events in the WIMP search data, for the defined 
nuclear recoil acceptance window. For each energy bin, 
the derived electron recoil rejection efficiency and the nu- 
clear recoil acceptance values are listed in Table Q] 

In addition to the statistical events leaking from the 
electron recoil band into the nuclear recoil acceptance 
window, we observed anomalous leakage events in the 



FIG. 2: Position distribution of events in the 4.5 to 26.9 keV 
nuclear recoil energy window, from the 58.6 live-days of 
WIMP-search data. (+) Events in the WIMP-signal region 
before the software cuts. (©) Events remaining in the WIMP- 
search region after the software cuts. The solid lines indicate 
the fiducial volume, corresponding to a mass of 5.4 kg. 

The 3D position sensitivity of the XENON10 detec- 
tor gives additional background suppression with fiducial 
volume cuts [2j|. Due to the high stopping power of LXe, 
the background rate in the central part of the detector 
is lower (0.6 events/keVee/kg/day) than that near the 
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Nuclear Recoil Equivalent Energy (keV) 



FIG. 3: Results from 58.6 live-days of WIMP-search in the 
5.4 kg LXe target. The WIMP search window was defined 
between the two vertical lines (4.5 to 26.9 keV nuclear recoil 
equivalent energy) and blue lines (about 50% nuclear recoil 
acceptance). 

edges (3 events/keVee/kg/day). For this analysis, the 
fiducial volume is chosen to be within 15 to 65 /is (about 
9.3 cm in Z, out of the total drift distance of 15 cm) 
drift time window and with a radius less than 8 cm (out 
of 10 cm) in XY, corresponding to a total mass of 5.4 kg 
(Fig. [HI • The cut in Z also removes many anomalous 
events due to the LXe around the bottom PMTs, where 
they happen more frequently compared to the top part 
of the detector. 

After all the cuts were finalized for the energy window 
of interest, we analyzed the 58.6 live-days of WIMP- 
search data. From a total of about 1800 events, ten 
events were observed in the WIMP search window after 
cuts (Fig. [3]) . We expect about seven statistical leakage 
events (see TablelJ) by assuming that the ALogio (52/51) 
distribution from electron recoils is purely Gaussian, 
an assumption which is statistically consistent with the 
available calibration data, except for a few "anomalous 
leakage events". However, the uncertainty of the esti- 
mated number of leakage events for each energy bin in 
the analysis of the WIMP search data is currently limited 
by available calibration statistics. Based on the analy- 
sis of multiple scatter events, no neutron induced recoil 
event is expected in the single scatter WIMP-search data 
set. To set conservative limits on WIMP-nucleon spin- 
independent cross section, we consider all ten observed 
events, with no background subtraction. Figure 2] shows 
the 90% C.L. upper limit on WIMP-nucleon cross sec- 
tions as a function of WIMP mass, calculated for a con- 
stant 19% C e ff, the standard assumptions for the galac- 
tic halo (2d| . and using the "maximum gap" method in 
[25| . For a WIMP mass of 100 GeV/c 2 , the upper limit 
is 8.8 x 10~ 44 cm 2 , a factor of 2.3 lower than the pre- 
viously best published limit [27| . For a WIMP mass of 
30 GeV/c 2 , the limit is 4.5 x 10 -44 cm 2 . Energy res- 
olution has been taken into account in the calculation. 




XENON10 (58.6 live days) 
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FIG. 4: Spin-independent WIMP-nucleon cross-section upper 
limits (90% C.L.) versus WIMP mass. Curves are shown for 
the previous best published limit (upper, blue) [27| and the 
current work (lower, red), assuming a constant 19% £ e //- 
The shaded area is for parameters in the constrained minimal 
supersymmetric models [rj. |29||. 

The largest systematic uncertainty is attributed to the 
limited knowledge of £ e // at low nuclear recoil energies. 
Our own measurements of this quantity (2lj | did not ex- 
tend below 10.8 keVr, yielding a value of (13.0 ±2.4)% at 
this energy. More recent measurements by Chepel et al. 
[22I ] have yielded a value of 34% at 5 keVr, with a large 
error. 

A comparison between the XENON10 neutron calibra- 
tion data and Monte Carlo simulated data, including the 
effects of detector resolution and energy dependence of 
C e ff, provides an effective constraint on the variation of 
£ e f f for all energies in the analysis range [28j . The con- 
stant C e f f assumption used to calculate the limits above 
shows reasonable agreement at the 10% level between 
the Monte Carlo predicted spectrum and the measured 
energy dependence and intensity of the single scatter nu- 
clear recoil spectrum. The C e ff assumption which gives 
the best agreement implies a slightly more sensitive ex- 
clusion limit, and is not quoted. A conservative exclu- 
sion limit was calculated by including estimates of possi- 
ble systematic uncertainty in the signal acceptance near 
threshold. Also included was an estimate of the uncer- 
tainty in the energy dependence of the neutron scattering 
cross sections used in the Monte Carlo simulations. The 
C e ff assumption which gives poorest sensitivity, while 
remaining consistent at the 1% level with the neutron 
calibration data, corresponds to exclusion limits as high 
as 10.4 x 10~ 44 cm 2 (5.2 x 10" 44 cm 2 ) for a WIMP mass 
of 100 GeV/c 2 (30 GeV/c 2 ). 

Although we treated all 10 events as WIMP candi- 
dates in calculating the limit, none of the events are likely 
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WIMP interactions. ALogio(S2/51) values for 5 events 
(compared with 7 predicted) are statistically consistent 
with the electron recoil band. These are labeled as No.'s 
3, 4, 5, 7, 9 in Fig. M and Fig. H As shown in Table [Q 
these leakage events are more likely to occur at higher 
energies. A posteriori inspection of event No. 1 shows 
that the ST coincidence requirement is met because of a 
noise glitch. Event No.'s 2, 6, 8, 10 are not favored as ev- 
idence for WIMPs for three main reasons. First, they are 
all clustered in the lower part of the fiducial volume (see 
Fig. [5]) where anomalous events happen more frequently, 
as discussed above. Second, the anomalous 5*1 hit pat- 
tern cut discussed earlier for the primary blind analysis 
was designed to be very conservative. An independent 
secondary blind analysis performed in parallel with the 
primary analysis, used a more stringent cut to identify 
anomalous hit patterns in 51 and rejected 3 (No.'s 6, 8, 
10) of these 4 candidate events. Third, the expected nu- 
clear recoil spectrum for both neutrons and WIMPs falls 



exponentially with energy, whereas the candidate events 
appear preferentially at higher energy. 

The new XENON10 upper limit on WIMP-nucleon 
spin-independent cross section further excludes some pa- 
rameter space in the minimal supersymmetric models Q 
and the constrained minimal supersymmetric models 
(CMSSM) (e.g. [HI). 
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